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Abstract 

The charge exchange reaction pd — > npp at 1 GeV projectile proton energy is 
studied in the multiple-scattering expansion technique. This reaction is considered 
{SJ . in a special kinematics, when the transfer momentum from the beam proton to fast 

%4 \ neutron is close to zero. The differential cross section and a set of polarization ob- 

servables are calculated. It was shown that contribution of the final state interaction 
. between two protons is very significant. 

in 
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1 Introduction 

> 

Within the several last decades the deuteron- proton charge exchange reaction has been 
studied both from experimental and theoretical point of view. A considerable interest to 
this reaction is connected, first of all, with an opportunity to extract some information 
about spin-dependent part of the elementary nucieon-nucleon charge exchange amplitudes. 
This idea was suggested by Pomeranchuk [1] already in 1951, but up to now it is actual. 
Later this supposition has been developed in [2, 3, 4]. It was shown, that in the plane- wave 
impulse approximation (PWIA) the differential cross section and tensor analyzing power 
T 2 o in the dp-charge exchange reaction are actually fully determined by the spin-dependent 
part of the elementary np — > pn amplitudes. The analogous result was obtained in [5], 
where this process has been studied in the Bethe-Salpeter formalism. 

The differential cross section of the dp — > npp reaction at 3.34 GeV/c deuteron beam 
has been measured in 1970's in the lm hydrogen bubble chamber of the JINR Syn- 
chrophasatron [6, 7]. However, the obtained statistics is not sufficient to evaluate the 
magnitude of the spin-dependent part of the elementary amplitude. Nowadays the exper- 
iment on the study of the dp-charge exchange reaction at the small transfer momentum 
in GeV-region is planned at ANKE setup at COSY [8]. The aim of this experiment is to 
provide the information about the spin-dependent np-elastic scattering amplitudes in the 
energy region, where the phase-shift analysis data are absent. 

From our point of view, under kinematical conditions proposed in this experiment, 
when momentum of the emitted neutron has the same direction and magnitude as the 
beam proton (in the deuteron rest frame), and relative momentum of two protons is very 
small, the final state interaction (FSI) effects play very important role. The contribution of 
the D-wave in the DWF into differential cross section in this kinematics must be negligible 
[9]. However, for the polarization observables the influence of the D-component can be 
significant. 

The goal of our paper is to study the importance of the D-wave and FSI effects under 
kinematical conditions of the planned experiment. We consider pd — > npp reaction in 
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the approach, which has been used by us to describe the pd breakup process at 1 GeV 
projectile proton energy [10]. This approach is based on the Alt-Grassberger-Sandhas 
formulation of the multiple-scattering theory for the three-nucleon system. The matrix 
inversion method [11, 12] has been applied to take account of the FSI contributions. 
Since unpolarized and polarized mode of the deuteron beam are supposed to be employed 
in the experiment, we also calculate both the differential cross section and a set of the 
polarization observables. It should be noted, in this paper we have not considered the 
Coulomb interaction in the (pp)-pair. This problem is nontrivial and requires a special 
investigation. 

The paper is organized as following. In sect. 2 the short description of the general 
theoretical formalism is given. The special kinematics, when the transfer momentum from 
the beam proton to neutron is close to zero, is considered in sect. 3. The results of our 
calculations for the differential cross section and polarization observables are presented in 
sect. 4. The figures in this section demonstrate the behaviour of this observables obtained 
in the PWIA and PWIA+FSI with D-wave in the DWF included and without it. The 
significant dependence of the calculation results on the elementary NN-amplitudes is also 
shown. We conclude with sect. 5. 

2 Theoretical Formalism 

In accordance to the three-body collision theory let us write the matrix element of the 
deuteron proton charge exchange reaction 

Pip) + rf (0) -> n(pi) + p(p 2 ) + p(p 3 ) (!) 

in the following form [10] 

U pd ^ npp = V2 < 123|[1-(2,3)][1 + t 23 (E-E 1 )g 23 (E-E 1 )}t^ 2 m mS)>, (2) 

where the operator g 23 (E — Ei) is a free propagator for the (23)-subsystem and the 
scattering operator t 23 (E — E\) satisfies the Lippmann-Schwinger (LS) equation with 
two-body force operator V 23 as driving term 

t 23 (E - E x ) = V 23 + V 23 g 23 {E - E 1 )t 23 (E - E 1 ). (3) 

Here E is the total energy of the three-nucleon system E = E\ + E 2 + E 3 . 

Let us rewrite the matrix element (2) indicating explicitly the particle quantum num- 
bers, 

Upd^npp = V2 < v\m X T X ,v 2 m 2 T 2 ,p 3 m 3 T 3 \\y - (2, i)\uj 23 t s ^ 2 m \pmT , ^iM d oo(23) >, 

where uj 23 = [l + t 23 (E — E\)g 23 (E — E\)] and the the spin and isospin projections denoted 
as m and r, respectively. The permutation operator for two nucleons was introduced 
here. The operator t^f™ is symmetrized NN-operator, t^f™ = [1 — (1, 2)]t 12 . Inserting the 
unity 

1 = / dp'IpWr' >< pm'r'l, 
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we get the following expression for the reaction amplitude 

J = (_i)i/2+^ < 1 -t 2 1 -t z \Tt 2 + r 3 X \t' 2 -t' 3 \Tt 2 + r 3 X \t\t' 2 \T 'M t , > 
< \t X - - t' 3 \T'M t , >< ^m 2 ^m 3 \SM s >< ^m' 2 ^m' 3 \SM' s > 



J dp > /p ,SM s 



t ST (E - Ei) 
l+m N — — 2 



Po',SM' s 



Po -Po +«0 

< p[ mi , (ph' + q/2) m' 2 \t T s ' ym {E - E' 3 )\pm, (p ' - q/2) m" > (4) 

< m"m' 3 \i> 1Md (po - q/2) > - (2 <-> 3), 

where E' 3 = \Jm? N + (p f — q/2) 2 , the nucleon mass , and we have introduced the 

momentum transfer q — p — pi, relative momenta p = \{p 2 — Pz) an d Po — \{.P~2 — Pz')- 
Henceforth, all summations over dummy discrete indices are implied. 

In momentum representation the DWF ipiM d (k) with spin projection M d is written as 

L 

NW£)>=£ E <LM L lM s \lM d >u L (k)Y^(k)\lM s >, (5) 

L=0,2 M L =-L 

with the spherical harmonics Yj^ L {k) and the Clebsh-Gordon coefficients in standard 
form. In our calculations, we have employed the following parameterizations of the S- 
and D- state wave functions 

u o(p) = \/-£ 2~T 2 > Mp) = \/-£ 32~ 2 ( 6 ) 

as proposed in refs. [13, 14, 15] . 

We assume that r = r 2 = r 3 = 1/2 and T\ = —1/2. Then the isotopic coefficient can 
be calculated and Eq.(4) is simplified 

J = \ < \™*\ m *\ SMs >< \™^m' 3 \SM' s > 

< LM L \M s \lM d >< ^m"^m' 3 \lM s > (7) 

/ dpo' < ^sm s tmM'SM' s TM t > uM - q/2\)Y^ L (fo ^2) 

< plrnx, (po' + g/2) m' 2 \t° sym (E - E' 3 ) - t\ ym {E - E' 3 )\pm, (p ' - q/2) m" > 
- (2<-3), 

The wave function of the final p-p pair 

< ^ { poSm s tm t \Po'SM' s TM t >= 8(p~o -Po')5 Ms m> + 

^ 2 "IV - n < PoSM s \t ST \p 'SM' s > (8) 
Po -Po +i0 

contains the FSI part, which can be taken in different ways. 

In this paper we use the matrix inversion method (MIM) suggested in refs. [11, 12] and 
applied to study the deuteron electro-disintegration [16, 17] and deuteron proton breakup 
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process [10]. As in ref. [16], we consider the truncated partial- wave expansion, 



< iPp SMsTM T \Po'SM' s TM T >= 5 Ms m' s $(po -Po) + 



(9) 



E E Y l ^p )<l P iSMs\JMj>^ ll (po')< l'n'SM l s \JMj>Y l 7\p '), 



J=0 Mj=-J 



where J max is the maximum value of the total angular momentum in p-p partial waves 
and a = {J, S,T} is the set of conserved quantum numbers. The radial functions ip^,(po') 
are related via 



to the partial wave functions (pf,i V (p' ) , which have the asymptotics of standing waves. The 
coefficients Ou» can be expressed in terms of the corresponding phase shifts and mixing 
parameters [16]. 

Within the MIM, the functions ipf v {p' ) can be represented as 



where the coefficients Bf[,(j) fulfill a set of linear algebraic equations approximately equiva- 
lent to the LS integral equation for the p—p scattering problem 1 . Here N is the dimension 
of this set, pj are the grid points associated with the Gaussian nodes over the interval 
[-1,1] and pn+i = Po (details can be found in ref. [18]). It should be noted that in this way 
the nucleon wave function is expressed by a series of ^-functions allowing one to reduce a 
triple integral in Eq. (7) to a double one. In addition, the method offers the opportunity 
to consider the nucleon wave function in the continuum directly in momentum space what 
simplifies all subsequent calculations. 

3 Collinear geometry 

In this paper we consider the special kinematics, when transfer momentum q—p — pi is 
close to zero. In other words, the neutron momentum has the same value and direction 
as the beam proton. In fact, since the difference between proton and neutron masses and 
deuteron binding energy take place, the transfer momentum is not exactly zero, q ~ 1.8 
MeV/c. But because of this value is very small and has no significant influence on the 
results, we shall suppose q = in the subsequent calculations. 

Under such kinematical conditions one can neglect the dependence of the high-energy 
nucleon-nucleon matrix Inn{E — E' 3 ) in Eq.(7) on the internal nucleon-nucleon momen- 
tum in the deuteron and express it in the center-of-mass system (c.m.s.) through three 
independent amplitudes 




(10) 




(11) 



?nn(q = 0) = A + (F - EXWX^Q*) + B{a 1 a 2 ), 



(12) 



We neglect here the Coulomb interaction in the (pp)-pair. 
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where 



\p* + p' | 

with the fast proton (neutron) momentum p* . We use some results of relativistic potential 
theory [19, 20] to relate this NN t-matrix in the c.m.s. with that in the frame of interest 
(see, also [10]). 

< mim' 2 , 1 1 \pp* ,mm" > = NN'F < m^D^u,^^ >< m' 2 \D\u,^ Q )\ii' 2 > (13) 

< Mi/4|*cm(p l ')|/V / >< n\D(u,p)\m >< [x"\D(u,jf )\m" >, 

where D is a Wigner rotation operator in the spin space and u is a four-velocity. In our 
kinematical situation, when p — p\ 3> f? , each of these operators is slightly different 
from the unit operator, so that in a good approximation the tATAr-matrix in the frame of 
interest has the same spin structure. The product of the normalization factors iV and N' 
and kinematical factor F is 

NN'F = (14) 

In such a way, we have following relation for high-energy NN t-matrix in the different 
frames of reference 



< 



mim' 2 |t(p,$,)|mm'' >= ^ ^ Ep < m 1 m' 2 \t cm (p*)\mm" > (15) 

To evaluate such quantities without their momentum angular decomposition we use 
the phenomenological model suggested by Love and Franey in refs. [21]. In this approach 
the corresponding matrix elements are expressed through the effective NN-interaction 
operators sandwiched between the initial and final plane-wave states, that enables us to 
extend this construction to the off-shell case. Obviously, such off-shell extrapolation does 
not change the general spin structure. 

Since the (pp)-pair belongs to the isotriplet, one can anticipate that the FSI in the x Sq 
state is prevalent at comparatively small p - va l ues - 111 suc h a wa Y we get the following 
expression for amplitude of the dp charge exchange process 

J = JpwiA + J^So 

Jpwia = mN + E ? < LM L IM S \IM D > u L (p' )Y L M ^pZ) 
'i hip 



{< 7,m' 2 ^m 3 \lM s >< m im2 \tl m (f) ~ tlm(f)\ mm 2 > 
1 1 

< ^rn' 2 -m 2 \lM s >< rnxm^t^Jj?) - t\ m {f ~)\mm' 2 >} 



(16) 



/ i \1— ni2- m' 2 _i_ rp 11 

J^so = { — 1 - A ^^<LM L lM s \lM D >5 m2 . m:i <-m >/ -~m > 2 \lM s > 

47T 4:hj p Z Z 

J dp 'p ' 2 J dp 'u L (p' )Y L M -(pi)^Z(Po) 
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One can see, we have very simple integral over angular variables p' . As the result of 
this integration we get following relation for Ji S{) 

(-l)i-^-"4 mjv + E 1 „1 

= 7^ 4£ m2 ~ ms 2 m 2 ~ m z\ 1MD > 



< mimi|^(T) " mm" > J dp 'po' 2 ^ZiP>oip' ) (17) 

Note, the integral over radial variable p' also has no any difficulties, since ^ooHK)) contains 
5-f unctions. 



4 Results and discussions 



We define unpolarized 2-^3 cross section by the standard manner 

p 



E r 

a (dp — > npp) = (27r) 4 — / dp\dp 2 dpz\ j\ 2 5 A (A — momentum), (18) 

P J 



where \J\ 2 = l/6Tr(J J + ) is the square of the process amplitude averaged over all 
particles spin states. Using the 5 3 -function to eliminate the p% integration and changing 
variables from pi to q we have 



a ( d p ^ npp ) = ( 27 r) 4 ^ [ dqdp 2 \J\ 2 b(m d + E p - Jm% + (p - q) 2 - E 2 - J m% + (q - 



E p 
p 

Taking p 2 ,q <C p this expression can be reduced to 

E 2 

a(dp -> npp) = (2tt) 6 ^| J dq 2 dp 2 dcosQ 2 p\\J\ 2 (20) 

We define general spin observable related with polarization of initial particles in terms of 
the Pauli 2x2 spin matrices a for the proton and a set of spin operators S for deuteron 
[22] as following 

_ Tr(Ja a S p J) 
Ca " ~ Tr(JJ+) ' (21) 

where indices a and (5 refer to the proton and deuteron polarization, respectively; cr 
and So corresponding to the unpolarized particles are the unit matrices of two and three 
dimensions. In such a way, Eqs.(16- 17) for dp- charge exchange amplitude enables us to 
get the relation for any variable of this process taking into account two slow protons final 
state interaction in -state. So, we have following expression for the spin- averaged 
squared amplitude 

Co ee Tr(^+) = ^(^^ (22) 

{F 2 - B 2 )w{p 2 ){w{p 2 ) - 2V2ReU(p 2 ))(3cos 2 9 2 - 1)}, 

where U{p 2 ) = u{p 2 ) + J dp' p 2 ijj ^ 1 (p' )u(p' ) lfi t ne S-component of the DWF corrected on 
the FSI of the (pp)-pair. One can see, the Cq is independent on the neutron and proton 
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azimuthal angles <p q and <p 2 . To obtain Eq.(20) we have considered this fact and performed 
the integration over azimuthal angles. We use a right-hand coordinate system defined in 
accordance to the Madison convention [23]. The quantization z-axis is along the beam 
proton momentum p. Since the direction of p x pi is undefined in the collinear geometry, 
we choose the y-axis normal to the beam momentum. Then third axis is x — y x z. 
The tensor analyzing power can be presented in the following form 

C , w -C = ^( !!! ^) 2 {4(F 2 -5 2 )(W 2 ( P2 ) +W 2 ( P2 )) + 

(2F 2 + B 2 )w{p 2 ){w{p 2 ) - 2V2ReU{p 2 ))(3cos 2 9 2 - 1) + (23) 
9B 2 w{p 2 ){w{p 2 ) - 2V2ReU{p 2 ))sin 2 6 2 cos2</) 2 - 
54(F 2 - B 2 )w 2 (p 2 )sin 2 9 2 cos 2 9 2 sin 2 <p 2 }. 

Note, that only squared nucleon- nucleon spin- flip amplitudes B 2 and F 2 are in expression 
for the tensor analyzing power C 0jyy and differential cross section. However, the spin cor- 
relation due to vector polarization of deuteron and beam proton contains the interference 
terms of this amplitudes 

C y ,yC = -^(^^) 2 i?e(W){2W 2 (p 2 )- W 2 ( P2 )- 

w(p 2 )(w(p 2 ) + V2ReU(p 2 ))(l - 3sin 2 6 2 sin 2 (j) 2 )} + (24) 
l2w(p 2 )Im(FB*)ImU(p 2 )(cos 2 9 2 - sin 2 6 2 cos 2 <p 2 ) 

It is interesting, that there is the term proportional to the imaginary part of U(p 2 ). It 
has a non-zero value only in case when FSI is taken into account. The analogous result 
we have obtained for the vector-tensor spin correlation 

C y , xz -C = -A(^^y/m(F5*){2W 2 (p 2 )- W 2 b 2 )- 

w{p 2 ){w{p 2 ) + V2ReU(p 2 ))(l - 3sin 2 6 2 sin 2 <j) 2 ) + (25) 
18w 2 (p 2 )sin 2 9 2 cos 2 9 2 cos 2 (f) 2 } — 
3V2w(p 2 )Re(FB*)l7nU(p 2 )( cos 9 2 — sin 9 2 cos (j> 2 ) 

In order to evaluate these observables we consider kinematics, when one of the slow 
protons is emitted along the beam direction as well as neutron, i.e. 9 2 = 0°. Under such 
conditions Eqs. (22)-(25) are significantly simplified. 

Co = ^(^^) 2 {(25 2 + F 2 )(W 2 ( P2 ) + ^ 2 (p 2 )) + 
(F 2 - B 2 )w(p 2 )(w(p 2 ) - 2V2ReU(p 2 ))}, 

Co, y yC = ^(^^) 2 {2(F 2 -fi 2 )(W 2 ( P2 )+ W 2 ( P2 ))+ (26) 
(2F 2 + B 2 )w(p 2 )(w(p 2 ) - 2V2ReU(p 2 ))} 
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Cy,yC = ~ ( mjV 2 ^ ) {Re(FB*)[2U 2 (p 2 )-2w 2 (p 2 )-V2ReL((p 2 )w(p 2 )]- 
3\/2Im(FB*)ImK(p 2 )w(p 2 )} 



C y , xz -C = -A ( mjV 2 ^ P ) 2 {Im(FB*)[2U 2 (p 2 ) - 2w 2 (p 2 ) - ^ifeMfoOwfo)] + 

3v / 2 J Re(F J B*)/mW(p 2 )w(p 2 )} 

The differential cross section and three polarization observables are presented in figs.(l- 
4). The Love and Franey parametrization with a set of parameters obtained by fitting of 
the modern phase shift data SPOO [25, 26] has been employed for NN-amplitude. The full 
lines correspond to calculations with taking into account both the FSI in the (pp)-pair 
and S- and D- waves in the deuteron. The results obtained in the PWIA are shown by 
the dashed lines. The dash-dotted lines in the figs. (2-4) for polarization observables are 
the full calculation results without D-wave in the DWF. In fig.l the full and dash-dotted 
lines are undistinguished. All calculations were carried out with Paris NN-potential [24] 
and Paris DWF [13]. 

One can see, the FSI contribution to the differential cross section (fig.l) is signifi- 
cant even at the very small proton momentum, while for the polarization observables 
the difference between PWIA and PWIA+FSI is visible only for p 2 > 10 - 15 MeV/c. 
However, with increase of the proton momentum up to 50 MeV/c the importance of the 
FSI corrections to the PWIA also increases. 

Note, the absolute value of the tensor analyzing power C§ m (fig. 2) in the momentum 
interval of interest is near zero. In order to understand the source of that, we disregard the 
D-wave in the DWF. Then, as a consequence from Eq.(26), the polarization observables 
are defined by the ratio of the nucleon-nucleon charge exchange amplitudes only 

I F 2 -B 2 
°> yy ~ 2 ' 2B 2 + F 2 

C - 2 Re{FB * ] (27) 



c —1 ■ 



Im(FB*) 
2B 2 + F 2 



Thus, the nearness of the tensor analyzing power to zero indicates that the absolute values 
of the spin-flip NN- amplitudes approximately equal each other, \B\ ~ 

The vector-tensor spin correlation C VjXZ (fig.4) has also very small value, \C y/xz \ rs 0.06. 
The magnitude of this observable decreases up to zero for p 2 50 MeV/c, if the FSI 
corrections and D-wave in the deuteron are taken into account, while it is almost constant 
in the PWIA and PWIA+FSI without D-wave. One can see from Eqs.(26, 27) for C y>xz , 
the reason of this behaviour is connected with the small value of the imaginary part of the 
nucleon-nucleon amplitudes product, Im(FB*). In such a way, the great contribution into 
C V)XZ gives the term proportional to Re(FB*), which defined by D-wave and imaginary 
part of the generalized function U(p 2 ). Note, that ImU(p 2 ) ^ 0, if FSI taken into account. 

The other situation is for the vector- vector spin correlation C y<u (fig. 3). The term 
proportional to Re(FB*) gives also a considerable contribution in this observable , but 
it is multiplied on the U 2 (p 2 ). The magnitude of C V:y is close to the theoretical limit 
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-2/3, that confirms to the conclusion about approximate equality of the nucleon-nucleon 
amplitudes, \B\ and \F\. Besides, this allows to conclude, that the relative phase between 
these amplitudes is close to zero. It is seen from Eq.(27), where D-wave was neglected. 

Since all the considered observables are defined by the elementary nucleon-nucleon 
amplitudes mostly, it is interesting to compare their behaviour for different NN- 
parametrization. In figs. (5-8) we present the same observables as in figs. (1-4) for two 
sets of parameters. The full line corresponds to the parameterization based on the mod- 
ern shift analysis SPOO [25, 26]. The dashed line is obtained using a set of parameters 
for NN-amplitude from [21]. One can see, the difference for the differential cross section 
(fig.5) is about 1.5-2 times. The absolute value for the tensor analyzing power C 0:yy (fig. 6) 
with new parametrization is about 2.5 times smaller than that with parametrization [21]. 
The opposite situation is for the vector- vector spin correlation C y>y (fig. 7), where new pre- 
dictions is 1.5 times lager in comparison with old parametrization of the NN-amplitude 
[21]. The predictions of these two parameterizations for the vector- tensor spin correla- 
tion C y>xz (fig. 8) are even in opposite sign. Nevertheless, the qualitative behaviour of the 
curves in figs. (6-8) for different sets of parameters is similar. 

5 Conclusion 

We have studied the deuteron -proton charge exchange reaction at 1 GeV energy in special 
kinematics, <f ~ 0. The influence of the D-wave in the deuteron and FSI between two slow 
protons has been considered. It was shown, that D-wave and FSI effects are negligible for 
the polarization observables at the proton momentum up to 10-15 MeV/c. As the result, 
in this region the polarization observables are defined by the ratio of the nucleon- nucleon 
charge exchange amplitudes only. However, it should not be ignored, that importance 
of the D-wave and , especially, FSI into polarization observables increases at p 2 > 15 
MeV/c. In such a way, we conclude, that ratio of the nucleon- nucleon charge exchange 
amplitudes and phase shift between their can be extracted from experimental data rather 
simple, if the experimental conditions and technical setup possibilities allow to work in 
this small momentum interval. In opposite case, this procedure is more complicated and 
model dependent. It should be remembered, that FSI contribution into the differential 
cross section is very significant in comparison with PWIA predictions even at very small 
proton momentum. This fact does not enables us to get the absolute value of the nucleon- 
nucleon spin flip amplitudes without considering the FSI corrections. 



We are grateful to V.V.Glagolev, M.S.Nioradze and A.Kacharava for inspiration of 
interest to this problem. The authors are thankful to V.P. Ladygin for fruitful discussions. 
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Figure captions 



Fig.l The differential cross section at q = as a function of one of the slow pro- 
ton momentum. The dashed and full line correspond to the PWIA and PWIA+FSI, 
respectively. 

Fig. 2 The tensor analyzing power C yy vs. p 2 - The dashed line corresponds to PWIA; 
dash-dotted and full lines are PWIA+FSI without D-component in the DWF and with 
it, respectively. 

Fig. 3 The spin-correlation C y<y due to the vector polarization of the deuteron. The 
curves are the same as in fig. 2. 

Fig. 4 The spin-correlation C ViXZ due to the tensor polarization of the deuteron. The 
curves are the same as in fig. 2. 

Fig. 5 The differential cross section at q = as a function of p 2 . The dashed and full 
line correspond to the full calculation with a set of NN-amplitude parameters taken from 
[21] and fit of SPOO data [26]. 

Fig. 6 The tensor analyzing power C yy vs. p 2 - The curves are the same as in fig. 5. 

Fig. 7 The spin-correlation C y , y due to the vector polarization of the deuteron. The 
curves are the same as in fig. 5. 

Fig. 8 The spin-correlation C y>xz due to the tensor polarization of the deuteron. The 
curves are the same as in fig. 5. 
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